Kondo effect with non collinear polarized leads: 
a numerical renormalization group analysis 
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The Kondo effect in quantum dots attached to ferromagnetic leads with general polarization 
directions is studied combining poor man scaling and Wilson's numerical renormalization group 
methods. We show that polarized electrodes will lead in general to a splitting of the Kondo resonance 
in the quantum dot density of states except for a small range of angles close to the antiparallel case. 
We also show that an external magnetic field is able to compensate this splitting and restore the 
unitary limit. Finally, we study the electronic transport through the device in various limiting cases. 

PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b, 73.23.Hk 



I. INTRODUCTION 

One of the most remarkable achievements of recent 
progress in nanoelectronics has been the observation of 
the Kondo effect in a single semiconductor quantum dot 
(QD).— i 2 - Since these early experiments, the Kondo ef- 
fect has also been observed in carbon nanotube QDs,^>£ 
and molecular transistors^ connected to normal metallic 
leads. The success of these experiments has opened the 
way to a systematic analysis, both theoretically and ex- 
perimentally, of the effects of the leads correlations on the 
Kondo effect. Of particular interest are the experimental 
realizations of the Kondo effect connected to supercon- 
ducting leads 7 - or to ferromagnetic electrodes.— The latter 
is particularly important since it is connected to the new 
emerging field called "spintronics",- One goal of spintron- 
ics is to provide devices where the transport properties 
are controlled by magnetic degrees of freedom. In this 
sense, a Kondo quantum dot connected to ferromagnetic 
electrodes may be regarded as an interesting setup to 
study interaction effects in spintronics devices. 

From the theoretical point of view, it has been pre- 
dicted, despite some early controversial results ) 10 : 11 ' 12 ' 13 
that ferromagnetic electrodes with parallel polarization 
produce a splitting of the Kondo zero bias anomal y 14 : 15 : 16 
whereas polarized electrodes with antiparallel polariza- 
tion generally do not i 14 : 16 : 17 Some of these predictions 
have been checked experimentally by contacting Ceo 
molecules with ferromagnetic nickel electrodes. 8 The un- 
derlying mechanism responsible of the splitting of the 
Kondo resonance is quite clear: virtual quantum charge 
fluctuations between the dot and the leads become spin 
dependent. As a result, an effective exchange magnetic 
field appears in the dot which generally lifts the spin 
degeneracy on it. For antiparallel configurations, the 
contributions of each electrode to the effective field may 
eventually cancel each other for symmetric couplings. 



Most of the previous works have focused quite surpris- 
ingly only on parallel and antiparallel polarizations. Ser- 
gueev et al.— considered non collinear polarizations but 
did not take into account the spin splitting of the dot en- 
ergy level which naturally gives rise to spurious results. 
In a recent preprint, Swirkowicz et a/.— analyze the non 
collinear situation using the equation of motion method 
supplemented by an effective exchange field that has to 
be determined self-consistently. Such method gives qual- 
itative results for temperature T ~ Tk , where Tk is the 
Kondo temperature. The main result is that the splitting 
of the zero bias Kondo anomaly decreases monotonically 
when increasing the angle 26 between the two lead mag- 
netizations. In this paper, we go beyond the equation 
of motion method by combining a poor man scaling ap- 
proach and the numerical renormalization group (NRG) 
approach. The poor man scaling provides us with an an- 
alytical expression for the Kondo temperature as a func- 
tion of the lead polarization p and 6 [see Eq. l|12p ] while 
NRG allows us to compute dot spectral densities and 
the near equilibrium transport properties like the con- 
ductance. The plan of the paper is the following: In Sec. 
Ull we describe our model Hamiltonian. In Sec. IIIII we 
use poor man scaling to estimate the Kondo tempera- 
ture. In Sec. IIV1 we present the NRG results for the dot 
density of states and for the conductance. Finally, Sec. 
fVl contains a brief summary of the results. 



II. MODEL HAMILTONIAN 



In order to describe this system we use the following 
general Hamiltonian 



H = Hl ea ds + H tun + Hdot, 



(1) 



2 



where Hi eads denotes the Hamiltonian of the conduction 
electrons in the leads and is defined by 



fe,7,± 



(2) 



In this equation ± creates an electron with energy 
£fc± in lead 7 = L, R with spin along ±n 7 where n 7 is 
a unitary vector parallel to the magnetization moment 
in lead 7. The angle between Hl and Hr is 29. Since 
the leads are assumed to be polarized, there is a strong 
spin asymmetry in the density of states Pl ±(w). The 
tunneling junctions between the leads and the dot may 
be described by a standard tunneling Hamiltonian 



(*7 C fe, 



(3) 



where d CT destroys an electron in the dot with spin a =f 
, J.. Here the possible values for the dot spin are quantized 
along the axis e z — (nj, + nn)/\n.L + nn\. The operator 



fc,7,<7 

and reads: 



c a ,fc,T 



is a linear combination of the operators c 



cos(%) -sin(%0 
sin(4) cos(^) 



Ca,k- 



(4) 



where Q = ±0 here. H dot is the usual Anderson Hamil- 
tonian and reads 



H, 



dot 



2J e d ,<jd\d s 



(5) 



where S z = (djd T - The last term of Ec l- © 

describes the Zeeman energy splitting of the dot. In 
what follows we mainly consider a symmetric situation 
in which t 7 = t. We denote the total tunneling strength 
by T = T L + T R with T L = nt 2 (p + +p-). The effective 
magnetic fields generated by the ferromagnetic electrodes 
is then along the axis e z = {nL+nii)/\nL+n]i\. We next 
perform a poor man scaling analysis. 



III. POOR MAN SCALING ANALYSIS 

The analysis extends in a straightforward manner the 
one developed in Ref. (T4|. In this part we neglect the 
energy dependence in the density of states and also sup- 
pose p 7 ± ~ p±. In the ferromagnetic leads, the ratio 
p = (p+ — p-)/ (p+ + p- ) is in general different from zero 
and as shown bellow it is one of the relevant parame- 
ters to describe the effect of the spin polarization. For 
an energy dependent density of states, the situation is 
more complex as will be shown in Sec. IIVI One per- 
forms a two stage renormalization group (RG) analy- 
sis by reducing the cutoff A from D the bandwidth to 
£d = £d] + £di which renormalizes the parameters of the 
Anderson model. In general, £j and ej, renormalize differ- 
ently due to the spin dependent density of states (DOS) 



m the leads M- The spin degeneracy of the Anderson en- 
ergy levels is therefore lifted with an effective Zeeman 
splitting 



5e a 



ed] - Sdi ~ pTcos(9) ln(D/A). 



(6) 



This splitting prevents from reaching of the strong cou- 
pling regime. Nonetheless, an external magnetic field B 
can be applied locally to compensate this internal field 
and restore the degeneracy between e-p and £j_. One 
can then perform a Schrieffer- Wolff transformation. The 
Kondo Hamiltonian so obtained is similar to the usual 
Kondo Hamiltonian: 



Hk = \ £ (J z 'c{ a1 c m - J zi c{ al c m )S z 

k.k',a,0 

+ J± ( c Li c wS+ + c\ aX c m S-). 
The bare values of the Kondo couplings are 



j 1 - = r A = j zl ~ 2t 2 



u 



(\e d \(e d + U) 



= J, 



(7) 



(«) 



where e d ~ £ d *\ ~ £d[- Note also that V'lt is a linear su- 
perposition of ipL+ and ipi,-. One can then perform the 
second stage of the RG analysis directly on the Kondo 
Hamiltonian. The RG equations can be obtained by in- 
tegrating out the high energy degrees of freedom between 
a cut-off Ao ~ e d and A and read: 



dX 



d\nl 



dX 



d\nl 
dX± 
dial 



cos 2 (6>/2) 
cos 2 (<9/2) 
cos 2 (9/2) + 
cos 2 (9/2) + 



2 



1-p 
1-p 



sm z (9/2) 
sm 2 (6>/2) 



Ai, (9) 



1+p 

sin 2 (0/2) 



1 ~P _ : . 2 



1+p 
1 +P.,. 2 



1-p 



sm\9/2) 



Ai, (10) 
X±X Z ^ 
X±X z i, (11) 



where we have introduced the dimensionless Kondo cou- 
plings A ± = ^Jp+p- J x , X z y l = p± J zA/l and I = A /A. 
The dimensionless Kondo couplings are all driven to 
strong coupling at the same energy scale Tk which can 
be determined explicitly by integrating out the above RG 
equations. The Kondo temperature can be expressed as: 



T K (p,9) = A exp 



arctanh(pcos(6>)) 
2Xop cos(9) 



(12) 



where Ao = Jo (p+ + P-Y This expression generalizes the 
one obtained in Ref. [14[ to all values of 9. The depen- 
dence of the Kondo temperature with p, 9 appears as a 
function of the single parameter pcos(9) and is therefore 
a uniform function of 9. Tk takes its maximum value T K 
independent of p at 9 = ir/2 which corresponds to an- 
tiparallel magnetizations and Tjc is minimum for 9 = 
as expected. 
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IV. NUMERICAL RENORMALIZATION 
GROUP ANALYSIS 

In order to compute dot spectral densities and trans- 
port properties we proceed with a NRG analysis. At this 
point it is convenient to perform two unitary transforma- 
tions on the initial Hamiltonian defined in Eq. ([I]). We 
first make an even/odd transformation introducing the 
operators 



(13) 



Ck,e/o,± = (Cfe,i,± ± C k ,R,±)/V2. 

We next introduce a new basis defined by 



h,e,f/l — Cfe, e ,± , bk,o,1/l 



Cfe, 



(14) 



The tunneling Hamiltonian can be written in the new 
basis as: 

H tan = (*A,<A + ff - c (15) 

fc,r/— e/o,o" 

where t e = ty/2coa(9/2), and t a = frv/2 sin(<9/2). The 
crucial point is that the angle 9 arising due to the non 
collinear lead magnetizations is now hidden in the tun- 
neling amplitudes t e and t . 

Note that in this formulation the dot is coupled with 
a single effective lead with an angle and spin dependent 
bulk spin DOS given by 



eff 



P± cos 2 (6»/2) +p T sin 2 (<9/2), 



(16) 



which is mixture of the initial bulk spin densities of 
states. This allows us to treat this problem in a sim- 
plified way using the numerical renormalization group 
metho d 21 ' 22 ' 23 with the generalizations to treat arbitrar- 
ily shaped densities of states.— 

We insist that such formulation has been made possi- 
ble only because we consider symmetric tunneling ampli- 
tudes tL = ta. Note that for antiparallel polarizations 
(9 = 7r/2) we have = p e ^ and the effective model 
has spin symmetry. 



A. Study of the dot density of states 

From hereon we will consider the DOS along the po- 
larization direction in each lead to be given by the tight- 
binding expressions 



P±(w) 




uj + Se± 
D 



(17) 



where 5e + = — <5e_ = Se is the spin-dependent shift of 
the bands with respect to the Fermi level. The quantity 
Se characterizes the magnetization P of the leads given 
by 



P 



duif(u})[p+(uj) - p-(uj)}, 



(18) 



where f{uf) is the Fermi function. Note that the densities 
of states in the leads are no longer constant which is more 
realistic to describe ferromagnetic leads.— For the sake 
of simplicity, from hereon we take the chemical potential 
equal to zero. Note that for this particular choice the 
parameter p defined in the previous section is zero while 
the spin polarization is not. For small Se and tempera- 
tures the magnetization is simply given by P ~ In 
the present case, the most important effect of the leads 
magnetism is the occurrence of an internal field B int that 
creates the effective Zeeman splitting already mentioned 
above. 

To calculate the dot density of states we need the local 
Green function. In the standard form, we have 



G- 1 M = G - 1 M-S(a;) ) 



(19) 



where S CT (uj) is the self-energy and G 0a (uj) 1 is the non- 
interacting Green function 



G <r(uj) = 



1 



w - £d : cr - gpBa - A a (uj) 



(20) 



The hybridization function A a (uj) is spin-dependent and 
has in general a non-zero real part that is given by 



K[A a (0)] = a4t 2 Secos(9)/D 2 



(21) 



for our choice of the DOS and uj ~ 0. This is equivalent 
to an angle-dependent effective magnetic field 



p B B mt =8t 2 Secos(9)/D 2 



(22) 



that polarizes the dot's spin and destroys the Kondo ef- 
fect for pBB m t > T K . 

The first thing we want to check is whether the polar- 
ized leads with a given angle 9 are able to split the dot 
DOS \p a (uj) = — X[G CT (w)]/7r]. It has been shown that po- 
larized leads with antiparallel polarized directions (cor- 
responding to 8 = n/2 here) do not split the dot density 
of states i 13 ' 14 Is there a finite interval around 9 = tt/2 
where this property holds? We have plotted in Fig. 1 
the dot spectral density p^ in the spin up sector for dif- 
ferent values of 9. As expected the low energy peak of 
Pi is pinned at u> = for 9 — but is shifted away from 
uj = as soon as 9 is decreased. This is much clearer in 
the right inset of Fig. 1 where the peak is zoomed. 

The splitting of the total dot density of states p = 
P\ + Pi is shown in the left inset of Fig. 1 for four dif- 
ferent values of 9. The splitting gradually increase as 6 
is decreased from 9 = n/2 and is maximum for 9 = 0. 
For very small deviations of 9 around tt/2, the effective 
magnetic field psBint < Tk and no splitting is expected. 

For collinear polarization directions, it has been shown 
that this splitting can be compensated by an external 
magnetic fielcU^ restoring the Kondo effect. For general 
directions of the lead polarizations, an effective magnetic 
field is also generated. This effective magnetic field can 
be compensated by applying an external magnetic field 
in the plane (rzL,^i?) with a direction opposite to e z . 
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FIG. 1: Dot density of states for the spin-up sector as function 
of energy for different values of 9 = 0, 0.4557T, 0.485-7T, n/2. 
The shift of the bands is fixed to Se — 0.5D. Other parameter 
are: U = 0.2, = — U/2, and t = 0.071. Right inset: zoom 
of pi for the same values of 9. We clearly see the shift of the 
maximum Pj when decreasing 9 away from 9 = n/2. Left 
inset: Total density of states as function of frequency for the 
four aforementioned values of 9. The density of states is split 
as soon as the effective magnetic field is larger than the Kondo 
temperature. 



We therefore expect that the splitting of the Kondo zero 
bias anomaly induced by the polarized electrodes can be 
corrected by an adequate magnetic field for all values of 
9. We have checked that this property holds using NRG. 
We have presented in Fig. [2] the DOS p\ as function of 
energy for 9 — (where the splitting is maximum) and for 
different value of the external magnetic field. For a given 
value of the external magnetic field, the splitting of the 
Kondo peak vanishes and the Kondo peak is restored at 
lo = 0. This is also true for a general angle 9 (see Fig. [21). 
These properties of the dot's spectral density will clearly 
reflect in the transport properties of the system as we 
will see. 



B. Transport properties 

Using the Keldysh formalism, it has been shown that 
for the collinear magnetizations, as in the case of non- 
magnetic leads, the finite temperature conductance can 
be put in terms on the dot's spectral density.— This is not 
the case for arbitrary angles. The problem is that for non 
collinear magnetizations, the two spin channel of the L 
and R leads are mixed and the standard approach used to 
eliminate the lesser propagator can not be applied. To 
our knowledge, in these circumstances, the only direct 
way of evaluating the conductance at any temperature is 
through the Kubo formalism. This requires evaluating 
two particle propagators, a matter that requires sophis- 
ticated numerical codes. There are however some limit- 
ing cases that illustrate the general behavior as we will 



FIG. 2: Dot density of states for the spin-up sector as function 
of energy for different values of the magnetic field in the case 
with parallel magnetizations (9 = 0). Other parameters as in 
Fig. 1. The spin polarization in the dot and the splitting of 
the Kondo peak can be compensated with a magnetic field. 
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FIG. 3: Same as Fig. [2] with 9 = 7T/ 



see. In our NRG formulation we have Pd-\{0) — Pdj,(0). 
Since the DOS in the leads are energy dependent (see Eq. 
(HU), this equality does not prevent to have a non zero 
polarization following Eq. lfl~8|l though p = 0. For this 
case, the hybridization matrices T aa > = — I(Y> aa >) (where 
(7, a' =|, |) are diagonal and angle independent. This al- 
lows the use of the Meir Wingreen formula.-^ The T = 
conductance is thus simply proportional to |Gd(0)| 2 . The 
T = conductance has been computed using NRG and 
is shown in Fig. [4] as a function of 9 for different values 
of the splitting Se. The conductance reaches the unitary 
limit for 9 — ir/2 (antiparallel case) and is minimum 
for parallel magnetizations. Though this result may look 
unusual, this is expected since the splitting is zero for 
9 = tt/2 restoring the Kondo effect and is maximum for 
9 = 0. Remember that p(0) = with the definition of 
the DOS we used in NRG. 

This result should be contrasted to the situation where 
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PIG. 4: Conductance as a function of the angle for differ- 
ent values of the bands shift 8e = 0.1D, 0.25D, 0.513. Other 
parameters as in Fig. 1. 



an external magnetic field able to restore the Kondo zero 
bias anomaly is added. We consider constant bulk DOS 
as in Sec. IIIII Electrons with spins | have a different 
DOS at the Fermi energy than electrons with spins |, 
resulting in p(0) ^ (as in [IHI1I!]). At T = 0, the 
Kondo singlet is therefore formed and the Fermi liquid 
theory can be applied. The conductance G through the 
device reaches at T = its maximum conductance G u 
which reads: 



2e 2 



1 -p 2 



h l-p 2 cos 2 ((9) ' 



(23) 



Notice that here G reaches the unitary limit for 9 = 
whatever the value of p and is minimum for 9 = ir/2 
at fixed p. In order to determine to determine G u , we 
have extended the scattering approach developed by Ng 
and Lee 2 ^ to polarized electrodes following [13]. This 



anomalous denominator comes from the fact the imagi- 
nary part of the dot retarded Green function G r dG is sim- 
ply — l/r CTCr at T = which depends on the angle 9 and 
p. At finite T -C Tx{p,9), deviations of order (T/Tr-) 2 
from G u are expected using the Fermi liquid theory. On 
the other hand, at high temperature T >• Tk{p,9), the 
conductance is much smaller and can be computed by 
renormalized perturbation theory: 



2e 2 



p 2 cos(20)) 



3tt 2 /16 
In 2 [T/T K ] 



(24) 



Notice that the expression for the high temperature con- 
ductance has the same angular expression as the one for 
a tunnel junction between two ferromagnetic leads con- 
versely to the low temperature case. 



V. CONCLUSION 

In this paper, we have studied a quantum dot in the 
Kondo regime between two non collinear ferromagnetic 
electrodes by combining poor man scaling and NRG tech- 
niques. We have shown using NRG that the Kondo zero 
bias anomaly is in general splitted for all angles between 
the two lead magnetizations (except for the antiparallel 
case). This splitting can be compensated by an external 
magnetic field. Based on the spectral densities, we have 
addressed transport properties for both compensated and 
uncompensated cases. 

Note added Some of the results presented in this pa- 
per have been also independently obtained in a recent 
preprint. 
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